Critical particle-hole composites at twice the Fermi wave vector in U(l) spin liquid 

with a Fermi surface 
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We find " chiral symmetry breaking" at finite energies in U(l) spin liquid, corresponding to critical 
particle-hole composite states with twice of the Fermi momentum (2fej?). We investigate this Fermi 
surface problem based on the Nambu-Eliashberg theory, where the off diagonal pairing self-energy 
is introduced to catch the Aslamasov-Larkin vertex correction. This approach is quite parallel with 
the case of superconductivity, where such Aslamasov-Larkin quantum corrections in the particle- 
particle channel are well known to be responsible for superconducting instability, formulated as the 
Nambu-Eliashberg theory in an elegant way. We obtain the pairing self-energy, which vanishes at 
zero energy but displays the same power law dependence for frequency as the normal Eliashberg self- 
energy. As a result, even the pairing self-energy correction does not modify the Eliashberg dynamics 
without the Nambu spinor representation, where thermodynamics is described by the typical z = 3 
scaling free energy. We discuss physical implication of the anomalous self-energy identical to the 
conventional Eliashberg normal self-energy, focusing on thermodynamics. 

PACS numbers: 71.10.-w, 71.10.Hf, 71.27.+a 



I. INTRODUCTION 

Non-Fcrmi liquid physics in thermodynamics and 
transport has been at the heart of condensed matter 
physics. Quantum criticality is regarded as one source of 
such emergent phenomenal, where exotica such as elec- 
tron fractionalization 2 and Landau-Ginzburg- Wilson for- 
bidden duality 3 beyond a microscopic model were pro- 
posed to explain unique but universal critical phenom- 
ena. 

An effective field theory for quantum critical dynam- 
ics can be generally written in terms of renormalized 
electrons interacting with order parameter fluctuations 4 , 
associated with either symmetry breaking or emergent 
gauge symmetry. Recently, it has been clarified that two 
dimensional Fermi surface problems are still strongly in- 
teracting even in the large- N (spin degeneracy) limit as 
the non-Abelian gauge theory^—, where one dimensional 
processes embedded in two dimensions, i.e., forward and 
backward scattering near the Fermi surface turn out to 
be relevant. One patch formulation was originally ex- 
amined to keep only the forward scattering, where Lan- 
dau damping for critical boson dynamics is shown to be 
exact 5 . Two patch construction was proposed to be a 
minimal model for the Fermi surface problem, taking 
into account not only the forward scattering but also the 
backscattering associated with 2kp particle-hole excita- 
tions, where kp is the Fermi momentum. The Landau 
damping dynamics turns out to be still preserved while 
the Eliashberg fermion self-energy gets corrected via so 
called Aslamasov-Larkin (AL) vertex corrections, pertur- 
batively performed to result in an anomalous exponent 
proportional to 1/N-. Critical boson dynamics associ- 
ated with SU(2) symmetry breaking (ferromagnetism) 
was considered, arguing that even the dynamical expo- 



nent may be modified 7 . 

Unfortunately, all previous calculations have been per- 
formed in a perturbative way based on the Eliashberg so- 
lution, thus it is difficult to discuss on non-perturbative 
critical dynamics of bosons and fermions. Recently, we 
performed an infinite order summation for particular ver- 
tex corrections, given by ladder diagrams^. This was 
performed in a fully self-consistent way, resorting to the 
Ward identity. It turns out that this class of quantum 
corrections does not modify the Eliashberg dynamics in 
the one patch formulation. 

AL vertex corrections were perturbatively demon- 
strated to play an important role in both boson and 
fermion critical dynamics, as mentioned before. In this 
respect we need to perform an infinite order summation 
for this class of quantum corrections. Physically speak- 
ing, these quantum processes are possible to cause non- 
trivial dynamics for 2kp particle-hole excitations. It is 
well known that the same class of quantum corrections is 
responsible for superconducting instability of the Fermi 
liquid state, where a pole with an imaginary frequency in 
the particle-particle t-matrix is identified with a Cooper- 
pair stateii. Similarly, 2kp particle- hole composites may 
result from AL quantum processes, expected to modify 
the critical dynamics. 

In this paper we propose how to introduce the AL 
vertex correction into the two dimensional U(l) gauge 
theory with a Fermi surface non-perturbatively. Resort- 
ing to the analogy with superconductivity, where the 
superconducting instability described by the AL ver- 
tex correction is reformulated by the anomalous self- 
energy in the Eliashberg framework of the Nambu spinor 
representation^, we argue that the off-diagonal self- 
energy associated with the 2kp particle-hole channel in- 
corporates the same class of quantum corrections in the 
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Nambu spinor representation, where the Nambu spinor 
is composed of two opposite Fermi momenta, kp and 
—kp. We evaluate an anomalous pairing self-energy 
in the Nambu-Eliashberg approximation, which vanishes 
at zero energy but displays the same power law de- 
pendence for frequency as the normal Eliashberg self- 
energy. As a result, even the pairing self-energy correc- 
tion does not modify the Eliashberg dynamics without 
the Nambu spinor representation, where thermodynam- 
ics is described by the typical z — 3 scaling free energy. 
However, nonzero pairing self-energy corrections at finite 
frequencies break " chiral symmetry" in the two patch 
formulation of the U(l) gauge theory, indicating criti- 
cal particle-hole composite states with 2kp. We discuss 
physical implication of the anomalous self-energy iden- 
tical to the conventional Eliashberg normal self-energy, 
focusing on thermodynamics. 

II. NAMBU-ELIASHBERG THEORY FOR THE 
ASLAMASOV-LARKIN VERTEX CORRECTION 



where f s „ represents a spinon field and a describes U(l) 
gauge field in the spin liquid state, s = ± represents the 
patch index, corresponding to + = +kp and — = —kp, 
respectively, and a = 1,...,N expresses the spin index. 
Note that the sign in the gauge coupling is different be- 
tween the two patches, implying the gauge-current min- 
imal coupling, rj is an infinitesimal coefficient to control 
artificial divergences in quantum corrections, which can 
be cured by self-energy corrections^, z is the dynamical 
exponent determining the dispersion relation of gauge 
fluctuations. It is given by z = 3 for several problems 
such as ferromagnetic or nematic quantum criticality in- 
cluding the present spin liquid problem 1 ^ while z = 2 in 
the spin density wave ordering 4 . Both the Fermi velocity 
vp and the curvature 1/m are set to one. Generally, one 
may see that this effective field theory describes dynamics 
of fermions interacting with collective bosons, associated 
with symmetry breaking. 



A. The ladder vertex correction 



We start from the two patch formulation for the two 
dimensional U(l) gauge theory with a Fermi surface 



J dr J d 2 r j /j^ ^?j<9 T — isd x — d^j f s „ Dynamics of bosons and fermions can be described self- 

o * ii consistently by their self-energy corrections. The fermion 

+ —=afl (T f s <j + a(—dy)^~a>, (1) self-energy is expressed as follows 
\/N y ) 



S s (A;o) = -i J J (!^2 A ( fco + Qo,k + q;ko 1 k)G s {k + q 0l k + q)D(q ,q) 



1 f dqo f d 2 q 



N J 2tt J (2tt) 



T^-lk F ( k o + qo,k + q; k , k)G- s {k + q 0l k + q)D 2 k F {qo,q), 



(2) 



where the first part results from the forward scattering respectively. H(qo,q) is the boson self-energy, basically 
associated with G s (ko + qa, k + q) while the second contri- given by the fermion polarization function 
bution comes from the backward interaction related with 
G_ s (fco+qo, ^+9): n °t introduced in the one patch formu- 
lation. G s (ko,k) is the fully renormalized Green's func- 
tion and D(q Q ,q) is the fully renormalized boson propa- 
gator, given by 

G s (k ,k) ' 



ir)k Q + sk x + fcg - T, s (k ) ' 



D(q 0l q) = 1 r— T 1 TT / v ( 3 ) 

\Qv\ + n (9o»9) 
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dko 



d 2 k 



. A(k a + qo,k + q;k ,k)G s (k + qo,k + q)G s (k ,k), 

Ztt j (Ztt) 



/dko f d 2 k 
~2~7 J Jr^ A 2k F ( k o + qo,k + q;k ,k)G s (ko + qo,k + q)G- s (ko,k), (4) 




where U2k F (qo,q) is the fermion polarization function 
around the 2fc_p transfer momentum. 



FIG. 1: Fermion self-energy [S s (fco)] and boson self-energies 
[H(qo,q) and Il2k F (qo,q)] '■ The fermion self-energy results 
from both forward and backward scattering and 2fci? polariza- 
tion should be taken into account explicitly in the boson self- 
energy. The thick line represents the fermion Green's function 
and the wavy line does the gauge propagator. The shaded re- 
gion corresponds to the ladder-type renormalized vertex. 



A(k + qo,k + q;k ,k) and A 2feF (fco + <?o, k + q; k , k) are 
vertex corrections, where the former is associated with 
the forward scattering and the latter is related with the 
backscattcring, given by 



A(fc + q ,k + q;k ,k) = 1 



J J (2tt) 2 + 1o-lo,k + q-l;k -lo,k- l)G s (ko -l ,k- l)G s (k + qo~h,k + q- l)D(l , I) 



(27T) 

■ dL r d 2 i 



~ J ~2n J (2w) 2 ^ ko + q ° ~ l ° ~ m °' k + q - I - m;k a - lo - mo, k-l-m) 
G s (k - lo - m ,k - I - m)G s (k + q - lo - m , k + q - I - m)D 2 k F (m ,m) 
G- S (k -lo,k- l)G s (k + q - la, k + q - l)D 2kF (l , 



(5) 



and 



A 2 fe F (fco + qo, k + q; k a ,k) = 1 



dl 



d 2 l 



. , . s A 2 k F (ko + q - lo,k + q- I; k - l , k - l)G- s (k - l , k - l)G s (k + q - lo, k + q - l)D(l ,l) 
Zir J (Zir) z 

dl 



A 2 fe F (fco + qo~lo- mo, k + q - I - to; k a - l - m Q , k-l-m) 

- m)D 2 k F (m ,m) 



2tt J (2tt) 2 

G- S (k -l - mo, k-l- m)G s (fc + q - lo - m , k + q 
G s (k -lo,k- l)G s (k + q - lo,k + q - l)D 2kp (l , 1), 



(6) 



r 



respectively, in the ladder approximation. It is straight- D(lo, I) to be G s (ko — lo, k — I) and scatters into G s (ko + 
forward to read these vertex corrections. The first part q a — l ,k + q+l) by A(k +qo — lo, k + q— I; k a — l a , k — l), ac- 
in the forward scattering vertex is that G s (ko,k) emits cepting D(lo, I) to be G s {ko+qo, k+q). In the second cor- 
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FIG. 2: (Color online) The ladder vertex correction near the 
zero momentum transfer turns out to be irrelevant in the 
Eliashberg solution. 




FIG. 3: (Color online) The ladder vertex correction near the 
2k f transfer momentum turns out to be singular in the Eliash- 
berg solution, and should be taken into account carefully. 

rection G s (ko, k) emits £*2fe F (lo, to be G- s (ko—lo, k—l), 
which scatters into G s (ko — Iq — rriQ,k — I — m) by 
D2k F ( TO o> m )- G s (k Q — l Q — m ,fc — I — to) becomes 
G s (k + q - l - mo, k + q - I - to) by A(fc + q - 
Iq — rrio, k + q — I — m; fco — Iq — mo, k — I — to), and it 
accepts D2k F (mo, m) to be G s (fco + qo — h, k + q — I). It 
again accepts Dik F (£oi I) to be G s (ko+qa, k+q). One can 
understand the 2kp vertex correction in the same way. 
Equations (2), (3), (4), (5), and (6) consist of a fully self- 
consistent framework for a particular class of quantum 
corrections in the Fermi surface problem. 

If the backward scattering given by A 2 fc F (fc + qo,k + 
q;ko,k) is neglected, these self-consistent equations are 
reduced to those of the one patch formulation. As men- 
tioned in the introduction, such forward scattering con- 
tributions do not modify the Eliashberg dynamics^. The 
ladder vertex correction was solved, resorting to the 
Ward identity with an ansatz, where the vertex function 
is given by the fully renormalized Green's function. Such 
an ansatz turns out to recover both the non-interacting 
limit and one dimensional case. Inserting this ansatz ex- 
pression into both self-energy equations, we obtain fully 
self-consistent solutions, basically the same as the Eliash- 
berg self-energies in the low energy limit. 

However, the role of the backscattering was empha- 
sized because the 2kp ladder vertex turns out to be 
singular—, given by 



Mk F (ko + q ,k + q;ko,k) 




where the transfer momentum q is already expanded near 
2fcp- e_F is the Fermi energy and c is a positive numerical 
constant of the order of one. Ao is the bare vertex, set to 
be 1 in this paper. It should be addressed that this singu- 
lar 2kp particle-hole ladder vertex was not obtained self- 
consistently. Instead, it was found from a perturbative 
summation up to an infinite order (parquet approxima- 



tion), based on the Eliashberg Green's function^. Un- 
fortunately, we do not know how to get a self-consistent 
vertex solution, where the Ward identity cannot be used 
for the 2k p channel. In this subsection we discuss several 
possibilities based on the previous perturbative vertex so- 
lution. 

■q is an anomalous exponent for the 2kF vertex func- 
tion, depending on the spin degeneracy and given by 
X/N basically. Inserting this vertex into Eq. (4), the 
previous analysis found two kinds of boson self-energy 
corrections^ 

y ujov f I Wo^ 2kF > 

/^Z{ c „(») 2 '-» + c lVpI|"-'}-', (8) 

y ujovp I Vw / 2kp > 

where the critical value is r\ c — 1/3. loq = (j^J , 

where vp and po are the Fermi velocity and curvature, 
respectively, and g is the gauge coupling constant. Re- 
member vp,po,g = 1 in this paper. 

The 2kp boson self-energy does not diverge in 77 < r\ c 
and can be neglected at low energies due to a positive 
constant V . On the other hand, it becomes singular 
when rj exceeds rj c , implying quantum criticality of the 
boson state. An important question is how the fermion 
self-energy behaves, scattered by this 2kp critical bo- 
son dynamics. Does it diverge to give rise to an in- 
stability to a density wave (charge or spin) or vanish 
with an anomalous exponent? The latter will be iden- 
tified with a quantum critical state, characterized by a 
diverging 2kp spin susceptibility. We believe that this 
question is not addressed appropriately, i.e., fully self- 
consistently even in the relativistic U(l) gauge theory, 
where the 1/N expansion is well defined. It was demon- 
strated that the ladder vertex correction results in singu- 
larity to the staggered spin susceptibility, based on the 
X/N expansion in the relativistic U(l) gauge theory^. 
Since the fermion self-energy vanishes with an anoma- 
lous critical exponent, the quantum criticality character- 
ized by the divergent staggered spin susceptibility can be 
regarded as a stable phase, sometimes identified with an 
algebraic spin liquid as one possible phase of the param- 
agnetic Mott insulator^. This statement is certainly plau- 
sible at N = N c , below which chiral symmetry breaking 
arises to result in an antiferromagnetic order, because 
jV = N c corresponds to an antiferromagnetic quantum 
critical points. Our real question is whether this critical 
solution survives or not away from the quantum critical 
point, i.e., in the case when N > N c . Although the pre- 
vious non-self-consistent analysis demonstrated quantum 
criticality of the antiferromagnetic spin susceptibility in 
the algebraic spin liquid state (N > N c ), it seems to be 
possible that a fully self-consistent treatment may cause 
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both the fermion self-energy and the staggered spin sus- 
ceptibility to vanish with non-trivial critical exponents. 

The present Fermi surface problem is much more com- 
plicated because the 1 /N expansion is not well defined 5 . 
Inserting the boson self-energy of Eq. (8) into Eq. (2), 
one can find the fermion self-energy, particulary corrected 
by the 2fcp process. As discussed above, it is not clear 
whether r\ exceeds its critical value or not. When rj is 
smaller than rj Cl the 2kp scattering is irrelevant due to 
gapped boson excitations at this momentum, and the 
Eliashberg dynamics will be preserved. However, the di- 
vergent boson self-energy in the case of r/ > r\ c makes 
the fermion self-energy more singular than the Eliash- 
berg one. Actually, inserting the singular 2fci? ladder 
vertex and the divergent 2fcp boson self-energy into the 
equation for the fermion self-energy, we obtain 

Wfco;^ c )oc-z^^M|fc |^, (9 ) 

which diverges in the k§ — ¥ limit. This can be inter- 
preted as a particle-hole bound state, analogous with a 
Cooper pair in the particle-particle channel, resulting in a 
density wave state, where gauge fluctuations play the role 
of a pairing glue, reflected in their divergent self-energy. 
We believe that this 2kp issue should be addressed 

I 



where g q is the electron-phonon coupling constant and 
r(po + Qo,P + IiPQiP) is the associated vertex. D(qo,q) 
is the phonon propagator. Here, the imaginary time for- 

I 



r(po + q ,p + q;po,p) = g q 

f dk f d 3 k . 



where T pp [(ko,k;p + qo,P + q){kQ + qo, k + q;po 7 p)} is the 
particle-particle t-matrix to show the transition process 
that (po,p) and (fco + qo, k + q) scatter into (po + qo,p + q) 



self-consistently in both relativistic and non-relativistic 
U(l) gauge theories. One way to resolve this issue is 
discussed in the last subsection. In this paper we assume 
that the 2k f ladder process does not cause any change 
to the Eliashberg dynamics, i.e., 77 < rj c , and investigate 
other quantum corrections. 



B. Review on the Aslamasov-Larkin vertex 
correction associated with the superconducting 
instability 

We review the role of AL quantum corrections in the 
Eliashberg dynamics, demonstrated to cause an anoma- 
lous exponent beyond the Eliashberg fermion self-energy. 
As discussed in the introduction, such AL vertex correc- 
tions in the particle-particle channel were shown to cause 
superconducting instability to the Fermi liquid stated. 
Interestingly, this t-matrix construction was reformulated 
in the particle-hole Nambu representation, resulting in 
the same critical temperature. In this subsection we focus 
on the connection between the Nambu-Eliashberg theory 
and the AL quantum corrected self-energy. 

The electron self-energy is given by 



+ q;Po,p)G(p + qo,P + q)g q D(q ,q), (10) 



I 

mulation is applied at zero temperature. The electron- 
phonon vertex can be approximated as follows (Fig. 4) 



qo,P + q){k + q ,k + q;po,p)}G(k + q Q ,k + q)G(k , fc), 

(11) 

I 

and (ko,k) with the transfer momentum (qo,q), respec- 
tively. Inserting this AL vertex into the fermion self- 
energy, we obtain 
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E(po) ~ - J J -^3 G (po + qo^p + q)g q D{q ,q)g q 

f dq f d 3 q 

G(p + qo,P + q)G(-p , -p)G(-p + q , -p + q)g q D(q ,q)g q , (12) 




where the most singular particle-particle channel was se- 
lected (Fig. 5). 




FIG. 5: (Color online) The electron self-energy by the AL 
vertex correction in the particle-particle channel for super- 
conducting instability 



The particle-particle t-matrix is given by the Bethe- 
Salpeter equation^ 



t pp[(-po,-p;po + qo,p + q)(-po + qo,-p + q;po,p)] = g q D(q ,q)g q 

/dkl f d 3 k' 
~2% J jr^ T pp^ Po + q °~ k 'o>~ p + q ~ k '> pa + k 'o> p + k ')(- p o + qo,-p + q^ 

G(po + k' >P + k')G(-p +qo~ k' , -p + q- k')g q - k ,D(q -k' ,q- k')g q - k >. 



See Fig. 4. Generally speaking, this integral equation 
cannot be solved analytically because it is not factor- 
ized. But, if the retardation effect can be neglected, one 
obtains its analytic expression^. 

We compare this approach with an elegant reformula- 
tion based on the Nambu-Eliashberg approximation. The 



(13) 

I 

Dyson equation is 

G(p ,p) = g(po,p) + g(po,p)^{po)G(p ,p), 
where the Green's function is a two by two matrix in the 
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Nambu spinor representation, given by 



G(p ,p) 



Cpt 



G(Po,p) F(Po,p) 
F*( Po ,p) -G(-po,-p) J ■ 



associated with particle-particle pairing. The electron 
self-energy is given by 



S(po) 



S(po) *(po) 

$*(p ) -s(-po) ; ' 



G(po 1 p) is a normal Green's function with a bare propa- where <&(po) is the pairing self-energy and T,(po), the nor- 
gator g(po,p) while F(p ,p) is an anomalous propagator mal one. Then, the Dyson equation becomes 



G{po,p) = g(po,p) + g(po,p)Z(k )G(p ,p) + g(p ,p)$(po)F*(p ,p), 
F*(po,p) = -g(-Po,-pW(po)G(p ,p) + g(-po,-p)Z(-p )F*(p ,p). 

I 



(14) 



Inserting 



F*(po,p) 



9{-Po,-p)®*(Po)G(p ,p) 
1 - g(~Po, -p)S(-po) 



from the second equation into the first, we obtain the full 
renormalized self-energy in the Nambu formalism 

y (n\-y(n\ a{-po,-p)Mvo)\ 2 

^tot(Po) - MPo] - t r, (i-o) 

i - g{-po, -p)S(-po) 

where the total self-energy is defined as 

G(p ,p) = g(Po,p) + g{Po,p)^tot(Po)G(p ,p)- 



The self-energy is determined in the Eliashberg approx- 
imation 



S (P°)=/ ^ / 7^?35pfc £) («)»9 , ) G! (P0 + «).P + 9)- 

(16) 



Inserting the Nambu Green's function into the above ex- 
pression, we find full self-consistent equations 



g(-po - go, -p - g)$*(po + go) 
- g(~Po - go, -p - q)^(-po - go) 



g(po + qo,p + q) 



i - g( P o + q o, P + q ){z(Po + go) - 1^™%^ } 



for the pairing self-energy and 



g(po + go, p + q) 



g( P o + q^ P + q){z( P o + go) - r^r-r/J^-ro) } 
I 



(17) 



(18) 



for the normal electron self-energy. imation and the Nambu-Eliashberg one, we perform the 

In order to see similarity between the i-matrix approx- following approximation for the normal self-energy 

I 



S(po) » f f -^j3 9(Po + qo,P + q)gphD(qo,q)g p h 

~ j " J J^f^( p ° + 9o ^ 2 ' 9 ^ + io,p + q)g{po + q ,p + q)g(-po - qo,-p- q)g P hD(q ,q)g ph . (19) 
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Comparing this expression with Eq. (12) in the i-matrix 
self-energy, we find some similarity if we identify the t- 



J 



matrix [Eq. (13)] with the square of the pairing self- 
energy [Eq. (17)] in the following approximation 



$*(po)«- 



dqo 
2ir 



d 3 q 
(2^ 



9lh D (Qo, q)g(po + qa,p + q)g(-Po -qo,-p- qWipo + qo)- 



(20) 



r 



It is clear that the i-matrix approximation is not equiv- 
alent to the Nambu-Eliashberg formulation in the normal 
state because the anomalous pairing self-energy should 
vanish as discussed in the next section. However, we 
believe that the anomalous self-energy in the Nambu- 
Eliashberg theory will catch an essential physics of the 
AL vertex-corrected electron self-energy at the critical 
point, where the singular part of the normal self-energy 
results from the pairing one in the Nambu-Eliashberg ap- 
proximation while from the i-matrix process in the AL 
vertex correction. One can regard the Nambu represen- 
tation as one possible way to select a better basis set for 
an instability channel via " unitary transformation" . 



C. Nambu-Eliashberg theory 

One can construct a self-consistent framework for bo- 
son and fermion self-energy corrections, introducing the 
AL vertex with the 2kp particle-hole t- matrix Bethe- 
Salpeter equation. However, this strategy is not practical 
because such self-consistent equations are much compli- 



J 



cated, where the retardation effect is inevitable. As the 
AL vertex correction is reformulated based on the Nambu 
spinor representation in superconductivity, we introduce 
such quantum processes in the same way, but the com- 
ponent of the Nambu spinor will be different from that 
in superconductivity. 

Resorting to the analogy with superconductivity, we 
rewrite the two patch model Eq. (1) in the Nambu spinor 
representation 

S e t! = [ dT j d2r { F ° fal - id x r 3 - dp) F a 
+^=aFlr z F a + a(-d^a}, (21) 



where F a 



f+(7 



is the Nambu spinor composed of 



each-patch fermion. x 3 acts on the +kp particle and — kp 
particle space. 

Both boson and fermion self-energies can be found in 
the Nambu-Eliashberg theory 



U(q ,q) =tr J 



dko 
~2T 

dqg 
2ir 



d 2 k 

(27T) 2 

d 2 q 



T 3 G(fc + 9o, k + q)T 3 G(k a ,k), 
2 T 3 G(k + q ,k + q)D(q , q), 



(22) 



r 



where 



G(k ,k) = (irjkol + k x r 3 + kp - S(fc )) 

is the fermion Nambu Green's function with the matrix 
self-energy 



$*(fc ) S(fco) 



S(fco) is the normal part and $(A)o) is an effective pair- 
ing potential. More explicitly, the normal fermion self- 
energy, the fermion pairing potential, and the boson self- 
energy are given by 
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S(fco) 
*(*b) 

n((7o,<?) 



1 f dq f d 2 q , , 

= ~Nj J^G(k + q ,k + q)D(q ,q), 

/dkn f fl 2 k ( 
2^ y ^{G(fco + 9o,fc + g )G(fc ,fc) + G(fco + go,-fc- g )G(fco,-fc) 



(23) 



where 



G(fc ,fc) 



iryfc - fc^ + k 2 - E(fc ) 



F(k ,h) = - 



(i V k + k 2 y - S(fco)) - fc2 - |$(fc ) 
*(fcb) 



(ir/fco + ^ - S(fc )) - - |$(fc ) 



(24) 



r 



arc normal and abnormal fermion propagators, respec- 
tively. 



D. A solution of the Eliashberg equation 

Although it is not easy to solve Eq. (23) in a gen- 
eral way, the Landau damping ansatz for the boson self- 
energy 



9o 



makes this problem straightforward. In appendix we will 
check that the Landau damping boson self-energy is self- 
consistent indeed. 



Inserting the Landau damping boson self-energy into 
the boson propagator, normal and anomalous fermion 
(25) self-energies satisfy 

I 



s<*°) = w 



dy 


y | 


/ 


dq 


7T 


y z + 1 J 


J-k a 


2tt 


dy 


V | 


f 


dq 


7T 


y* + 1 J 


J — ko 


2tt 



E(fc 



qo ) 



V / S2(fco+ < 7o)-|$(fco+9o)| 2 ' 
$(fco + go) 

v /s 2 (fc + < zo)-|$(fc + go)| 2 ' 



(26) 



These self-consistent equations are typical in the Nambu- 
Eliashberg theory. 

Considering the symmetry E(fc ) < — > <&(fco) in Eq. 
(26), we propose 

A 2 

E(fe ) = -i— sgn(fc )|fco|*, 

$(fco) = -^ S gn(fco)|fc |* (27) 



as a possible solution. Inserting this expression into Eq. 
(26), we see 
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Equating this with Eq. (27), we obtain 

(28) 

where A and \ cannot be determined, but only A 2 + x 2 
is given by Eq. (28). This is the trace of the S(fc ) < — > 
$(fc ) symmetry in Eq. (26). 

Appearance of the off diagonal self-energy is far from 
triviality. If the self-consistent equation for the anoma- 
lous part is linearized, one obtains the matrix equation 

A ko = M kok ' o A k ' o , 

where A ko corresponds to the pairing self-energy and the 
column index is frequency from zero to cutoff. This ho- 
mogeneous equation does not have a nontrivial solution 
generally except the case when the matrix M ko k' satis- 
fies the particular condition of Det(Af — pi) = 0, where 
p is the eigen value set to be 1. Such a condition can 
be satisfied only at T < T c in the case of superconduc- 
tivity, and the corresponding eigen vector A ka vanishes 
identically when T > T c . The present situation is quite 
analogous with the T = T c case, i.e., the critical point, 



I 

where Ak ^o — > but Ak ^o ^ 0. The off diagonal self- 
energy correction may be regarded as one characteristic 
feature of the spin liquid state with a Fermi surface. 

An interesting result from the anomalous self-energy 
is that the chiral symmetry U = exp(i#T3) is broken at 
nonzero frequencies, indicating critical particle-hole com- 
posites with 2kp . Figure 6 shows an imaginary part of 
the anomalous propagator in Eq. (24). As frequency 
increases, it is enhanced to a certain frequency and de- 
creased, following the self-energy behavior. The maxi- 
mum point becomes larger and its width gets broadened 
as the momentum k x of fermions increases. Since the 
anomalous pairing self-energy is identical to the normal 
self-energy, the normal Green's function displays quali- 
tatively the same feature as the anomalous propagator. 
Appearance of critical particle- hole composites with 2kp 
is a novel feature of the two dimensional U(l) gauge the- 
ory with a Fermi surface. 



E. Thermodynamics 

We find physical implication of the present study from 
thermodynamics. Thermodynamics can be investigated 
in the Luttinger-Ward functional approach^ 



/72 7 
JL^rjlnjTsG-^feo, k)) - V(k )T 3 G(k , k)} 

+ T E / ^{^(D-\ qo ,q)) ~U( qo , q )D(q ,q)} 

+ T E T E / J J^^tr{T 3 G(ko + qo,k + q)D(q ,q)T 3 G(k ,k)}, (29) 



where the last term is called the Luttinger-Ward func- Using the Nambu-Eliashberg equation, one can sim- 

tional, resulting from the Nambu-Eliashberg approxima- plify this expression and obtain the following free energy 
tion in this paper. 
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1 0.6 




FIG. 6: (Color online) An imaginary part of the anomalous 
propagator in Eq. (24). The thick (red) line corresponds to 
k x =0.1 and k y = 0.0 while the dashed (blue) line does to 
k x = 0.3 and k y = 0.0. The Landau damping coefficient is 
7 = 0.1. Increasing the damping coefficient 7 from 0.1 to 0.5, 
the hump becomes broadened slightly as shown in the inset 
figure. 



F(T} . - T £ £ £ £ * ln {(,; + 4 8g „ (4o)|4ol ») 2 - 4 - £ w »} 



A2 
~2tt 



2tt 



,M 
T^wT 



9 » + 7 ra 



■}• 



(30) 



r 



where A is a momentum cutoff. It is trivial to observe from the simple Eliashberg approximation, given by 
that this free energy is qualitatively the same as that 

I 



r° dky f 00 dk x , r /,, 



feo 

00 dg : 



r 



v + ^ 



7 k„l 



A B(fco)|fco|^) 2 -^} 



(31) 



I 9k I 



except the anomalous self-energy part. In this respect 
thermodynamics of the Nambu-Eliashberg theory is es- 
sentially the same as that of the Eliashberg theory. 

To perform integrals in Eq. (30) is straightforward. 
The boson sector gives rise to the conventional T 2 con- 
tribution while the fermion part results in the typical 
scaling expression F(T) oc T^ d+Z ^ z with z — 3 and d = 2 
in the low temperature limi t 17 ' 18 . 

An essential point in this demonstration is that ther- 
modynamics is determined by the z — 3 quantum criti- 
cally. This doe not seem to be consistent with the pre- 
vious perturbative analysis 6 although this study did not 
examine the thermodynamic potential in their perturba- 
tive scheme, thus a direct comparison with the present 
result is not possible. To construct the Luttinger-Ward 



functional consistent with this perturbative scheme is an- 
other problem. But, we expect that the singular part of 
the thermodynamic potential will be given by 



F a (T) cx 



-y 

k 



d 2 k 



ln[ 



-Gf{iko,k)]: 



where Gf(iko, k) is the full fermion propagator with the 
AL-type vertex corrected self-energy. Since the fermion 
self-energy has an anomalous exponent, the specific heat 
coefficient will follow scaling different from z = 3. 

One may claim that this difference is natural because 
the Luttinger-Ward functional Eq. (30) does not seem 
to include the AL-type vertex corrections. However, this 
is not true. The key point in the present analysis is 
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that the AL-type vertex corrections are not only intro- 
duced but also summed up to an infinite order in the 
Nambu-Eliashberg framework, taking into account the 
anomalous pairing self-energy. As a result, the anoma- 
lous exponent in the fermion self-energy disappears and 
thermodynamics remains qualitatively the same as that 
of the simple Eliashberg theory because the pairing self- 
energy is essentially the same as the normal self-energy 
consistent with the z = 3 scaling. 

It should be pointed out that although we argued that 
the singular contribution for such vertex corrections will 
be kept within this anomalous pairing self-energy at least , 
we could not prove the equivalence exactly and showed it 
approximately as discussed in section II-B. If this "equiv- 
alence" turns out to be true, our statement has impor- 
tant physical meaning. In this respect it will be valuable 
to perform a direct summation for the AL-type vertex 
corrections, expected to confirm our result. 

III. DISCUSSION AND PERSPECTIVES 
A. Discussion and summary 

Our study was motivated from a recent perturbative 
analysis^, based on the Eliashberg theory as the zeroth 
order, where the Aslamasov-Larkin (AL) -type vertex 
corrections give rise to an anomalous novel exponent for 
the fermion self-energy while they do not modify the 
Landau-damping boson self-energy. If this is a general 
feature beyond this level of approximation, various quan- 
tum critical phenomenal should be reconsidered because 
the emergence of the anomalous exponent is expected to 
result in various novel critical exponents for thermody- 
namics, transport, and etc. In this respect it was crucial 
to investigate the existence of the anomalous exponent 
beyond this level of approximation. 

Unfortunately, it was not easy to perform a direct non- 
perturbative summation for the AL-type vertex correc- 
tions due to retardation effects of gauge fluctuations al- 
though one can construct self-consistent equations analo- 
gous to superconducting instability. Instead, hinted from 
superconductivity (section II-B), we reformulated the 
non-perturbative summation, introducing an anomalous 
self-energy in the Nambu-spinor formalism, where the 
two components are momentum kp and — kp fermions, 
respectively. As a result, we found that the anomalous 
exponent of the normal self-energy turns out to disappear 
and the pairing self-energy has essentially the same func- 
tional form as the conventional Eliashberg self-energy for 
frequency. We argued that this is the result from the 
non-perturbative summation for the AL-type vertex cor- 



rections in a " unitary" transformed basis (Nambu-spinor 
formalism) . 

One implication of this result could be found in ther- 
modynamics. If one evaluates the thermodynamic poten- 
tial in a perturbative way with the AL-type vertex cor- 
rections, the appearance of the anomalous exponent in 
the fermion self-energy may give rise to a novel exponent 
for temperature dependence in the specific heat coeffi- 
cient, different from the z — 3 scaling result. Here, z is 
the dynamic critical exponent. However, we argued that 
the non-perturbative summation for the AL-type vertex 
corrections, performed in the Nambu-Eliashberg frame- 
work, leaves the specific heat coefficient unchanged from 
the conventional Eliashberg theory with z = 3. This is 
important for physical application to the heavy fermion 
quantum critical point, particularly, for the Kondo break- 
down quantum critical point scenario 1 ?, where the z = 3 
quantum criticality gives a qualitatively good descrip- 
tion for thermodynamics in a certain heavy fermion 
material^. 

More direct implication of the anomalous self-energy 
may be found in spectroscopy measurements. For exam- 
ple, if the 2kp susceptibility (spin or charge) is measured 
(or evaluated), we expect that the anomalous propaga- 
tor contributes and this will show unexpected composite 
states from critical fermions (kp and — kp) directly. Its 
detailed shape, for example, the peak or hump-like posi- 
tion and the width (life time of the state), is not clear at 
present. An optical conductivity (in the case of nematic 
quantum criticality) will also show an interesting feature. 
As expected, it would be enhanced at finite frequencies 
due to breaking of such preformed particle-hole pairs. 

We also discussed the role of 2k p ladder vertex correc- 
tions in the critical state described by the simple Eliash- 
berg theory. In particular, the issue was to prove exis- 
tence of such a critical state as shows power law diver- 
gence of the 2kF spin susceptibility while the fermion self- 
energy vanishes with a nontrivial exponent. Although 
this physics is typical at quantum criticality, we argued 
that this issue should be addressed more sincerely in the 
critical spin liquid state, i.e., self-consistently even in the 
relativistic U(l) gauge theory, where the 1/N expansion 
is well defined. 



B. Perspectives 

Although the Nambu-Eliashberg theory does not give 
an anomalous exponent beyond the Eliashberg approxi- 
mation, one can extend this framework, introducing the 
ladder-type vertex correction 
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/dko f d 2 k 
~2tt J Jr^Mko + qo^ + q-tko^^Giko + qo^ + qjrsGiko^), 

1 f dqo f d 2 q 

S(fco) = -— J — J —^A(ko + qo,k + q;k ,k)T 3 G(ko + qQ,k + q)D(q ,q), 



r 



(32) 



where the two by two matrix vertex function is given by 
the ladder approximation 



J 



/dil f d^l 
2tt J (2tt)2 A<yk ° 

T 3 G CT (fc + q - Iq, k + q - l)T 3 G a (k -l ,k- l)D(l , 1). 

I 



+ q a - l a , k + q - I; k -l ,k- I) 



(33) 



This can be regarded as natural generalization of the one 
patch formulation, taking the Nambu representation to 
incorporate 2k p backscattering. An important point is 
that the vertex function satisfies the Ward identity auto- 
matically because the Nambu-Eliashberg approximation 
is conserving, i.e., constructed from the Luttinger-Ward 
functional approach. We expect that the same strat- 
egy as the one patch formulation can be applied to the 
Nambu-Eliashberg theory for the two patch construction, 
where the Ward identity will help us introduce such ver- 
tex corrections^. Another advantage in this formalism 
is that the 2kp ladder vertex will be introduced natu- 
rally. It is an interesting future direction to solve these 
equations. 



We would like to thank T. Takimoto and S.-S. Lee for 
helpful discussions. K.-S. Kim was supported by the Na- 
tional Research Foundation of Korea (NRF) grant funded 
by the Korea government (MEST) (No. 2010-0074542). 



Appendix A: Evaluation of the polarization function 

We prove that the Landau damping boson self-energy 
Eq. (25) is preserved in the Nambu-Eliashberg theory. 
Inserting both normal and abnormal fermion Green's 
functions into the equation of the boson self-energy, we 
obtain 



J 



!%),<?) = 2 



dk 



-90 ^ n •> ■ 



dk 



v_ 

2ir 



2tt 



k x {k x + q x ) + k 2 y {k v + q y Y - i±\k \-k 2 y + i±\k a + g |T(fc y + q y ) 2 + *-j$~\ko + g | « |fcp| « 

{(k x + q x ) 2 - {ky + q y Y - 2zA|fc + qQ m ky + qy )2 + ^d|fc + qQ \i }{ fc 2 _ fc 4 + 2% X N § fc 2 + ^ ^ § y 



(Al) 



r 



where both normal and pairing self-energies are explicitly 
used. 



It is straightforward to perform the integration of k x , 
resulting in 



J 



n(g ,g) 



2?r 



dko 
~2^ 



dk„ 



Vd({B + C){C -D) + (B- C)q 2 x ) + Vc((B + D)(C - D) - (B - D)q 2 x 



2tt 



C^D({C-D) 2 -2{C + D)ql 



qi 



(A2) 
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where 



A2 A 2 A 2 + y 2 22 

B = kl(k y + q v ) 2 -l-\ko\hl + l-\k + qo\ i (k y + q y ) 2 + ^ |fc + go | 1 \ko\ 1 , 

A 2 A 2 + x 2 4 
C = (k y + q y ) 4 + 2i — \k + q \~ (k y + q y ) 2 \k + q \* , 

D = k y -2i^\k \hl-^^\k \i. (A3) 



Scaling k y and k as follows 

fc y = QyV, fc o = got, 



the above expression reads 



(A4) 



tt/ v 1 |?o| dt 



dy 

2 ?r \ly\ J-i 2lT J -00 2tt 



+ c)(c - d) + (b - c) J) + ^((6 + d)( c - d) - (6 - d) J) 
^Vd((c-d) 2 -2( C + d)| + |) 



(A5) 



where 6, c, and d are given by 

A 



b = q t{v 2 (y + 1) 2 - 4^ v K¥tT + l P + + 1)2 ( 



\ly\ 



IgolM 2 , A2 +X 2 



| (+1 ,i| t ,i(Wf)*} S4 ; 6 



AT2 



o - + 1/ + 24 |t + 1|% + lf (MI) 2 _ + 1|4 (Mi)'} . 9jc , 



Igolh 2 A 2 + X 2 ul i/M^ 



AT 2 



i<i*(W* 



r 



(A6) 



An important point is that this integral is convergent we obtain 
near y w 0. As a result, the main contribution comes 
from 



\v\ » i- 

Scaling y in the following way 



Y 



(A7) 



c - d w 4F J 



6 - d w 2y 3 

(A9) 



in the dominant region for the integral. Finally, we find 
(A8) the boson self-energy 



n(g ,g) 



1 M [° dt 



dY 



y 2 {4!--2y3|(Mf)- 4 } 



koj * 

dY 1 



2tt 



tt }_ x 2tt 



r 4 {i6y 6 -4F 4 |(Mf)-%|(Mf)- 6 } 



(A10) 



nothing but the Landau damping form. In this respect tion of the Nambu-Eliashberg equations. 
Eq. (25) and Eq. (27) are the fully self-consistent solu- 
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